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Abstract. It is suggested that the accretion disk excites the eigen modes of Alfven oscillations of the magnetic field 
tubes the ends of which are frozen to the neutron star surface. The resonance takes place when the eigen Alfven 
frequency coincides with the difference of the Keplerian disk rotation frequency and the frequency of rotation of a 
neutron star. This model explains the two narrow lines of QPO oscillations in kHz band and permits to calculate 
the masses and radii of the neutron stars in X-ray binaries. 
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1. Introduction 

Rossi X-Ray Timing Explorer (RXTE) in 1996 discovered 
the quasi-periodic oscillations (QPO) of X-ray emissions 
in low-mass X-ray binaries (LMXB) i n the band of 1 kH z 
frequencies ( van der Klis et al. 1996| ; van der Klis 2000| ). 
, They are narrow, the ratio of the line's width to the cen- 
^ tral frequency can be as low as 10^^, and in most cases 
5^ , there were observed two lines. The existence of such lines 
means that the accretion disk approaches very close to the 
surface of the rotating neutron star. Indeed, the Keplerian 
frequency 

where M is the neutron star's mass and r is the distance 
from the star's center, is of the order of 1 kHz at the 
distance of 1.7 • 10^ cm for the typical neutron star mass 
M ~ 1.4 • Mq. But estimated neutron star radius R is 
of about 10^ cm, and we see that the inner edge of the 
accretion disk can approach almost to the star surface. If 
we estimate the kinetic energy of a proton near the star's 
radius (r ~ R) it turns out of 3 • 10~^ erg or 20 MeV. For 
the characteristic proton's density in the accretion flow 
near the star surface, corresponding to the accretion rate 
M ~ 10" g/sec {L^ ~ ^Edd), rip ~ 8 • lO" cm^^^ we 
obtain the energy density of the accreting plasma of 

£ ~ 2.5 • 10^^ erg/cm^ 
Send offprint requests to: Ya. N. Istomin 



It seams that for the existence of the disk in the region 
near the star's surface, the star's magnetic field should not 
exceed the value of 



<£, or B < 2.5 • 10^ G. 



But this low value of magnetic field on the neutron star 
surface contradicts with the magnitude of the magnetic 
field observed for the milliseconds pulsars B ~ 10^~^° G, 
which are considered originate from the LMXB. But these 
estimations are incorrect. The plasma destroys the mag- 
netic field structure when its thermal energy exceeds the 
magnetic field energy, but not when its kinetic flow en- 
ergy exceeds the magnetic one. Moving across the mag- 
netic fleld lines plasma polarizes so that the electric fleld 
across the magnetic one appears. In the case of a thin disk 
the radial electric fleld 



c 



(2) 



is enough to carry the flow across the magnetic field with 
the Keplerian velocity — 2'KVY,r. For the dipolc mag- 
netic field {B (X r~^) the disk of the temperature T de- 
stroys the star's magnetic field at the distances 



r > rA = R 



B'^R^ 



1 1/4 



(3) 



4M(T/TOp)i/2 

If we substitute to (||) , the characteristic energy of X-ray 
emission ('^ 1 keV) as the maximum disk temperature we 
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Fig. 1. The accreting disk approaches so close to the neu- 
tron star surface, that plasma flow along the magnetic 
field lines. 



obtain 



r > 20R 



Vl09 G 



1/2 



Thus, we see that in the region where QPO of 1 kHz fre- 
quency can originate the structure of the magnetic field of 
the star is not disturbed by the accreting plasma. We can 
consider that the structure of the magnetic field near the 
star surface(r < lOi?) is given by the interior of the neu- 
tron star and has the dipole configuration. Besides, the 
accretion changes the initial neutron star magnetic field 
B ~ 10^2 G up to B ~ 10^ G. And it is naturally to con- 
sider the direction of the dipole magnetic momentum /x 
coinciding with the axis of star's rotation Jl. 

We assume that the accreting disk approaches close to 
the neutron star surface. Plasma fiows from the disk to 
the neutron star along magnetic field lines and form so- 
called "magnetic field plasma tubes" , which are shaped 
as magnetic field lines, are filled with plasma and have a 
small thickness. 

The narrow structure of the QPO of 1 kHz frequency 
suggests the resonance nature of its origin. We have the 
idea that it is a resonance between Keplerian rotation of 
the disk and an eigen modes of the oscillation of magnetic 
field tubes, frozen to the conducting star surface. 

Such model is more close to the sonic point beat fre- 
quency model (Miller et al. 1998), where the resonance 



takes place not with the Keplerian rotation of the disk, 
but with the radial motion of accreting fiow coinciding 
with the sonic velocity. The other model based on the 
eigen modes of the disk in the rotating frame was pro- 



posed by Titarchuk & Osherovich (1999). There is more 
exotic mechanism proposed by Klein et al. (199^ ) where 
QPOs are the oscillations of photon's bubbles. 



2. Magnetic tubes oscillations 

In this section we discuss the oscillations of the magnetic 
field lines with plasma frozen to them. Discovering this 
model it is convenient to use the dipole coordinates {I, a) 
introduced by the definition: 



B 
B 



[Va, V( 



(4) 



I is the coordinate along the magnetic field line 

r3 ^(f?2+z2)3/2' 

a numerates the magnetic field surfaces 



a = fL 



Q 



(g2_^ ^2)3/2 



const is the equation for the magnetic field lines 



^2 ^ (^)2/3^4/3 

a 



Q 



/" • 2, 
r — — sm ( 

a 



(5) 



(6) 



(7) 



Here (g, z) are the cylindrical coordinates, (r, 9) are the 
spherical ones. The value of a is also connected with the 
maximum deflection Qmax of the dipole magnetic field line 
from the star's center a = ^/ Qmax- For different magnetic 
surfaces the coordinate a changes from to ^/R. The co- 
ordinate / changes for a given a, from — (/i^ — iJ.RaY^'^ /R"^ 
to {ii^ — ^Ra)^/^ / R^ . The point / = corresponds to the 
top of the magnetic held line. Three coordinates (/, a, (j)) 
{(j) is the azimuthal angle) are the curvilinear orthogonal 
coordinate system. 

Because in the near (r < lOi?) region B^/Stt ^ £ 
the plasma can be considered as cold, the plasma density 
is proportional to the magnetic held strength, — XB. 
Here A is the constant not depending on the coordinates. 
The equations of the ideal MHD in rotating with angular 
velocity frame {fl is the star rotation frequency) can be 
presented for the electric current j as follows 



n ) 1 

rot rot <^ — 
LU^B I B 



B[jB] 



(8) 



where va is the Alfven velocity 

v\ = B'^/iTTUpmp = B/AnXrUp. 

We consider here that all the perturbed values j, SB, 
E depend on the time t and the azimuthal angle (j) as 
exp {—iLut + im(j)} . 

The electric current can be presented through its co- 
variant components ji, j2, j3 

J - Ji V / + j2 V a + j3 V (j), 

for which from (8) we obtain 



J2 1 - 



up- 



d 



uj^g^B dl 



dl 



d 



J3 



v^B 



Ji = 



u^B 



92 
dd 



2 O3S) 



Lu^g^B da 
d 1 d 



im 



dl dl 
v\g^B d 



{hB) = 

U2B) 



da 



9 9 r^9 d , . ^. im d , . ^. 

o-/B-QlinB)^-,^^{j,B) 



(9) 

(10) 

(11) 
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The boundary conditions for |l^, ^T]) arc on the star 
surface r — R . Due to the ideal conductivity of neutron 
star matter the tangential components of the electric field 
have to be zero. The electric field of the disturbance is 
proportional to the electric current 



iB 



{h V a + js V ( 



It gives the boundary conditions for j2 and 

32 = j3 = at r = i?, or 

I = ±lm = ± {^^^ - t^Raf^ / . (12) 

Equation describes the Alfven type oscillations be- 
cause it contains only the second derivatives over longitu- 
dinal coordinate I , but equation ( |lO| ) is that for the mag- 
netosonic waves (also the derivatives over a coordinate). 
In the case of nonaxial symmetric oscillations (m ^ 0) 
these modes are coupled. We consider the case when the 
accretion disk excites the MHD waves in the magneto- 
sphere with the dipole magnetic field. Rotating with the 
Keplerian frequency cjk different from the star rotation 
frequency 17, the conducting disk deflects the magnetic 
field lines in equatorial plane to 0-direction. Because this 
deflection is symmetric over (p we consider here only the 
case m = 0. As Alfven and magnetosonic oscillations are 
independent and Alfven wave has only the component 
of the disturbed magnetic field 



J2 



dSBs 



Ang^ dl ' 
SB = 6Bi V 1 + 5B2V a + 5B3V < 



SB:, 



g ■ SBa,, 



while the magnetosonic has two, SBi and SB2 
cB^g^ /dSB2 d5B{ 



J3 = 



47r V dl 



da 



(13) 



For this reason we choose only the Alfven wave {SBi 
SB2 = 0, SB^ ^ 0) 



J2 



d 



0. 



(14) 



uj^g'^B dl 
The boundary conditions 

j2 (a, I = ±^m) = 



imply the quantization of eigen frequencies of the Alfven 
modes 

uj = ujj{a), j = l,2,... 

After parameterization the eigen frequencies can be pre- 
sented in the form 



where 



(15) 



AnR^Xr 



1/2 




1.5 



2.5 



/ 



Fig. 2. The dependence of Alfven first harmonics on 
/• 



is the characteristic Alfven frequency, and tOj is the dimen- 
sionless amplitude depending on the dimensionless param- 
eter / = gmax/^, which means the height of the dipole 
magnetic field line in the star radius units. The calcula- 
tion for the general mode j = 1 gives the dependence of 
on / as pictured on Fig. 0. At large / the eigen frequency 
falls down like luj^ oc f^^'^- The reason of this behavior is 
that the main contribution to the frequency gives the wave 
structure in the top of magnetic field tube. Indeed, multi- 
plying the ( p^ ) on to (^B^g'^uj^/v\) • ^2 and integrating it 
over I from— to +lrm we get 



f-g^\B§iBj2)\'dl 



(16) 



From ([T^) it is seen that uj ~ ("a),^;^ /L, where 
L is the length of the magnetic field line. Because of 
VA = {B /ATrXmp)^^"^ oc f~^^'^ and L cx / we obtain the 
calculated dependence 

uj^r^'^ (/»i). 

As for the values of / near the star surface |/— 1| <C 1, then 
L (X (f — ly/^ and lu oc {f — 1)"-^/^ (/ ~ 1) as imagined 
on the Fig. ||. For the same reason the more high modes 
have the frequencies very close to the harmonics of the 
general mode 

ujj~j-uji, WA=t^i- (17) 

For the plasma density of rtp ~ 8 • 10^^ cm^"^ near the 
star surface and magnetic momentum ~ 10^^ G-cm"^ 
{B ~ 10^ G) the frequency wa is of the order of 27r kHz at 
the distance r ~ 5.6i?. It means that QPO of 1 kHz can 
appear as the result of the resonance of the disk rotation 
and Alfven eigen modes excited by the disk 

i ■ [uK - ri) = UJj ~ j ■ UJA- (18) 

Here i is the integer number. The resonance condition ( [l8| ) 
defines the magnetic field tubes, where the resonance oc- 
cur, : 



J 



(19) 



4 



Ya. N. Istomin and A. V. Makarov: Resonance accretion 




Fig. 3. There are two general solutions /i and /2 of the 

The solution of ( [l8| ) is pictured on the Fig. ^. In general 
case we have two solutions /i and /2 (/2 > /i)- Their 
positions are not constants in time because of the changing 
of the accretion rate and as a consequence of changing of 
the plasma density Up . The frequency of the Alfven eigen 



modes i 



-1/2 



. It means that higher accretion rate (or 



rip) larger /i and smaller /2. Accordingly wk(/i) becomes 
larger, but uj^if^) becomes smaller. The observed changes 
of QPO frequency correspond to the behavior of the root 
/i . Apart, it is observed the clear proportionality between 
intensity of X-ray emission and the frequency of QPO of 
1 kHz band. We consider that the resonance at / = /i is 
more strong than at / = /2 because of 

— it is more close to the star surface where the dipole 
magnetic field is large; and 

— usually /2 is large enough, so that cl'k(/2) falls down 
from the 1 kHz frequencies. 

Further we will discuss only the high frequency so- 
lution / = /i . How the QPO frequencies are connected 
with discussed resonance? At first the excited Alfven oscil- 
lations will destroy the accretion disk. The magnetic field 
gives the drift velocity Sv^ for the disk plasma particles 



5v, 



[Er SB^] 



SVz 1= VK 



SB, 



(20) 



(see Er from (H)). The disk matter will be pushed from 
equatorial plane in the direction along the dipole magnetic 
line to the star surface. And the accretion flow will be 
modulated by the frequency ^j{fi) + 0,. The point f = fi 
will be factually the inner edge of the accretion disk. In 
observed X-ray emission we can see several characteristic 
frequencies: 

1. ajK(/i), it is the modulation of the X-ray fiux by the 
inhomogeneity of the inner edge of the disk; 

2. (jjKifi) — i^, it is the beating frequency corresponding 
to the scattering of X-rays from the sources on the star 
surface, rotating with Q , by the inner edge of the disk; 

3. ijJj(fi) + il = i • Wk(/i) + (1 — j) ' ^1 it is the frequency 
of the accretion flow modulated by Alfven oscillations; 




/ 



1.5 2 2.5 3 

Fig. 4. The dependence of sonic first harmonics lo'^ on / 

4. different combinations of frequencies 1, 2, 3. 



3. Sonic oscillation 

It is naturally to consider that the plasma flow along the 
tube f — fi is not stationary. The sound mode can be 
excited. 

Let us consider plasma as isothermic with temperature 
T. For disturbance of concentration 6n along the magnetic 
tube we receive from hydrodynamics 



T 



■ divgradJn = 0. 



Sn 



In dipole coordinates it looks like 



T 



Q2 



(21) 



(22) 



The boundary condition is on the star surface. The plasma 
velocity component which is normal to the surface must 
be equal to zero. 



T 



V = — i- 



V Sn 



0. 



(23) 



i=±i 



The first harmonics corresponding to the solution of 
this equation with boundary condition can be presented 



^(/) 



^ UJs 



(24) 



where uj^ = T/nipR^ and uj^if) is dimensionless function 
on /, which is presented on the Fig. ^, Lui.{f) is of the order 
of 1 when / '--^ 2. In the region of / « 1 the function cu'sif) 
can be approximated by the expression 2/ (tt-^// — 1). So, 
ujQ is the charachteristic frequency for sonic oscillation. If 
we take T = 1 keV and R = 10 km then ujq = 30 
or fo = 4.9 Hz. Thus, sonic oscillation can be found on 
10 Hz band that is observed simaltaneously with the QPO 
of 1 kHz band. 
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Fig. 5. Due to Kerr metrics the disk precesses with Lense- 
Thirring frequency. 

4. Lense-Thirring precession 

While the accreting disk is spinning at the distance r 
about the neutron star with Keplerian frequency, it pre- 
cesses about its equihbrium positio n with Lense-Thirring 
frequency due to Ke rr metrics ( Thorne et al. 1986 ; 
Morsink fc Stella 19991) 



+ {a/cf (rg + r) ' 



(25) 



where a — J /M ( J — in is the star's angular momentum 
and / — 2ME? /7 is the momentum if inertia of the star), 
Tg — 2GAI/c^ is the gravitational radius of the neutron 
star of mass M and radius R. 

In the system of units where M is measured in masses 
of sun and R in 10 km the Lense-Thirring frequency can 
be presented in the form of 



WLT = 8.44 • 10" 



R 



(26) 



For M/R ~ 1, / ~ 1 and ^2 ~ 27r • 300 s^^ uj is of 
the order of 27r • 25 s^^. The second har monics of Lense- 
Thirr ing precession is usually observed ( ^tella fc Victri 



1999). In our model it can be illustrated by the fact that 



some mass is thrown out the disk each time it is in extreme 
position, i.e. with frequency 2ajLT (see Fig. 

5. Scorpius X-1 
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Fig. 6. Each line on this graph corresponds to one mea- 
surment of sonic and Keplerian frequencies. T is measured 
in eV, M in Mq, i? in 10 km. 
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We ha'je applied our speculations for Sco X-1 ( van der Klis 



et al. 1996 ). There were observed 4 QPOs: 6-20 Hz nor- 
mal/flaring quasi-periodic oscillations; 45 Hz QPOs and 
two kHz QPOs. According to our model we ascribed the 
6-20 Hz QPOs to the sound oscillations, 45 Hz QPOs to 
Lense-Thirring oscillations, and two kHz QPOs to Alfven 
and Keplerian oscillations correspondingly. 

Below we use the system of units where temperature 
is measured in eV, masses of objects - in masses of sun 
and distances - in 10 km. Thus, we can write 



Fig. 7. The variations in kHz QPO p eak separation as a 
function o f the lower kHz frequency ( Mendez fc van dei 
Klis 1999| ). It is well-seen that for Sco X-1 the dependence 



for sonic: a;,(/) = 0.98 \/ u;^(/). 



(27) 



for Keplerian: uo^if) - 11.5 • 10^ \ ^ /"^/^ (28) 

V R'^ 

where / = r/R - a parameter that characterizes the dis- 
tance between the disk and neutron star surface. 



is linear. 



For every pair of data (c^^k, ^^s) we can build a curve in 
M T \ 

— ; — ) space. This results in the graph on Fig. 0. 
i?'^ R J 

The axis of abscissae is the space for possible mean- 
ings of the ratio of M/R^ . For each of possible meanings 
of this ratio we can find the temperature T (or, to be more 
accurate, the ratio T/R?) for each experiment on deter- 
mining the frequencies and Wg. So we can evaluate the 
range of the change of temperature T . One can find it re- 
markable that all the curves cross each other in almost one 
point! So we can suppose that the temperature of the disk 
is approximately constant. From the graph we determine 
M/i?3 = 0.374 ± 0.004 and T/i?^ = 16.8 ± 0.8. 
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Fig. 8. Diagram of the dependence of gravitational mass 
on the radius for different equations of state of nuclear 
matter (Baym & Pethick 1979) and our point for Sco X-1 
on it. Larger neutron star radii are given by stiff (TI, MF) 
rather than soft (R, tt) equations of state. The classifica- 
tion corresponds to|Alcock et al. (1986). 



Mendez & van der Klis (1999) have observed that the 
separation between the two high-frequency QPO peaks 
changes from one observation to another. In our work we 
considered the equation 



i ■ (WK - = LOj ~ j ■ UJA- 

In general case the peak separation will be 



(29) 



n-ujA = { --i]i^A + ^- (30) 



It is easy to see, that when i ^ j the peak separation 
will not remain constant. For example, if i = 2 and j = 1 
then wk — = ^'-^a/2 + fi. If i = 1 and j — 2, then 
wk — = wa + Here we must mention that Alfven 
eigen frequency depends on the parameter A that changes 
when / will change and, therefore, when / changes, num- 
bers i and j can change. As a result, the dependence of 
the peak separation on ojk will be a piece linear function. 
The dependencies of separation frequency on the lower one 
for different objects is shown on the Fig. 0. For Sco X-1 
the dependence can be well approximated by the linear 
function (|^ with j/i = 3/4 and = 27r • 450 s'^ 

So, comparing the 45 Hz QPO with Lense-Thirring 
precession we determine the ratio M / R k, 0.60. Finally 
M = O.75M0, R = 12.6 km, T = 27 eV. The point (M, R) 
on the gravitational mass - radius diagram is shown on 
the Fig. I 

For this M,R,T we can find that wqs — 4.16 s^^ or 
vos = 0.66 Hz. This low value can be explained by the 
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Fig. 9. The X-ray color-color diagram for Sco X-1 
( |van dor Klis et al. 1996| ). Soft color is the (3-5 keV)/(l- 
3 keV), hard color the (7-20 keV)/(5-7 keV), count rate 
ratio. The bold line corresponds to the black body radia- 
tion with temperature varying from 1.7 keV to 1.95 keV. 



fact that disk approaches too close to the neutron star 
(down to / = 1.01) and in this region of / function ui'sif ) 
takes values much greater than 1 (e.g. Wg(l.Ol) = 15.7). 

At first sight, T must be equal to the radiation temper- 
ature of the neutron star, which can be determined from 
the X-ray color-color diagram. We put it on one graph 
from Scorpius X-1 and from black body (see Fig. P). The 
normal branch of that from Scorpius almost coincide with 
that of black body. And the temperature of the black 
body varies from 1.7 keV up to 1.95 keV. But if we es- 
timate the free path of the photons from the neutron star 
(l) ~ l/aup ~ 20 km we find that it is greater then R and 
so, radiation does not heat the inner part of the disk. 
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